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Abstract 

We study numerically the Sherrington-Kirkpatrick model as function of the magnetic 
field h, with fixed temperature T = 0.6T C . We investigate the finite size scaling behavior 
of several quantities, such as the spin glass susceptibility, looking for numerical evidences 
of the transition on the De Almeida Thouless line. We find strong corrections to scaling 
which make difficult to locate the transition point. This shows, in a simple case, the 
extreme difficulties of spin glass simulations in non-zero magnetic field. 

Next, we study various sum rules (consequences of stochastic stability) involving 
overlaps between three and four replicas, which appear to be numerically well satisfied, 
and in a non-trivial way. 

Finally, we present data on P(q) for a large lattice size (N = 3200) at low temper- 
ature T = 0.4T C , where, for the first time, the shape predicted by the RSB solution of 
the model for non-zero magnetic field is visible. 
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1 Introduction 



The Sherrington-Kirkpatrick (SK) model was introduced some time ago pQ as a mean field 
model for spin glasses. Its analytical, replica symmetry breaking (RSB), solution [2| displays 
in the glassy phase intriguing features such as an infinite number of pure states, described by 
an order parameter which is the non-trivial probability distribution of the overlap between 
two of them, P(q). 

The applicability of the mean field picture to short range spin glasses [3] is however still 
debated, and an alternative (family of) scenario called the droplet model has been put forward 
by several authors. One may in principle distinguish between these two theories of finite 
dimensional spin glasses by looking at the fate of the glassy phase for non-zero magnetic field. 
In the SK model (and accordingly in the mean field picture) one finds jH |S] that a magnetic 
field (of absolute value) lower than the critical De Almeida Thouless (AT) value 1iat(T) 
does not destroy the spin glass ordering, since the number of states is still infinite. On the 
other hand, in the droplet picture one has only two states at h = 0, related by the global 
inversion symmetry of the spins, and any small magnetic field makes the system paramagnetic 

mm 

It turned out unfortunately JU]-[2I] that strong finite size corrections make difficult to 
obtain a clear answer from equilibrium simulations in magnetic field (of the system sizes 
that can currently be handled). Very recently, for example, the authors of a study of the 
local excitations of the Edwards-Anderson model (EA) [23] claim that there is no transition 
whereas results on the out-of-equilibrium behavior of this model j21] appear in good agreement 
with the mean field picture. On the theoretical side, it has been recently proposed [22] that 
the transition below the upper critical dimension d u = 8 is governed by a fixed point different 
from the AT mean-field fixed point. 

This state of affairs motivated us to revisit the case of the SK model in magnetic field. 
There have been indeed few numerical studies of the SK model in magnetic field at fixed h 
[TBI Hp] . [2H]-j2n]; and the behavior of the system as a function of the magnetic field at fixed 
T has never been numerically investigated to our knowledge. 

We present results on the spin glass susceptibility and different dimensionless ratios of 
P(q) moments, looking for the quantities that are most appropriate for obtaining numerical 
evidence of the transition on the AT line. We will show in particular that finite size effects 
are strongly reduced if one consider the probability distribution of the absolute value of the 
overlap, and not the overlap itself. 

We moreover consider the overlaps between three and four replicas, checking the validity of 
some relations (the so called stochastic stability sum rules) which are an evident manifestation 
of the non-self-averageness of P(q) [301 131] that have been recently derived under very general 
properties 
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Finally we present data for P(q) in magnetic field, where the shape predicted by the 
solution of the model, with two peaks separated by a continuum is visible, for the first time. 
All simulations presented up to now show only a broad peak around q EA - Both large system 
sizes (we have N = 3200), and low temperature (we go down to T = 0.4) are needed to see 
this asymptotic shape. 



2 Model and Observables 

The Sherrington-Kirkpatrick spin glass model with TV sites jH E] is described by the Hamil- 
tonian 

Hj= ■ , '.i c7 ' n j ~ h J2 °i» (!) 

l<i<j<N l<i<N 

here cij = ±1 are Ising spins. The sum runs over all spin pairs and Jjj are quenched identically 
distributed independent random variables with mean value = and variance 1/N. We 
take Jij = ±N' 1 / 2 . 

In order to sample the probability distribution of the overlap P(q), one usually considers 
two independent replicas {a{\ and {t{\ evolving contemporaneously and independently: 

1 N 

JV i=l 



P(q) = Pj(q) = (5(q-Q)), (3) 

where the thermal average (•) corresponds to the time average during the simulation, whereas 
(•) stands for the average over the realizations. This is the order parameter of the model, 
which in the thermodynamic limit behaves like 

8(q- qEA) \h\ > h AT {T) 

P(q) = { x m 5(q-q m ) + P(q)+x M 5(q-q EA ) < \h\ < h AT (T) (4) 

\ \P(q) + P(-q)] + x ~f [5(q - q EA ) + 5(q + q EA )\ h = 0,T<T c 



where Ii A t(T) is the critical value of the magnetic field on the AT line, with Iiat{T) ~ 
(4/3) 1 / 2 (T c - Tf' 2 for T -> T 6 7 (T c = 1 in this model) 0. For T -> T~ one finds that 
x m oc q m oc /i 2 / 3 , and (q EA — q m ) oc (xm — %m) °t (h>AT(T) — h). Note that at h — the 
function P(q) is symmetric, reflecting the global flip {<jj} — > {— Uj} symmetry of the system, 
and the 5-function at q rn disappears. As we are going to discuss in detail, this implies a 
singular behavior for different quantities in the h — > limit, which is among the main sources 
of difficulties in finding evidence for the phase transition in magnetic field. 
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In order to locate numerically phase transitions, in our case the AT line, it is common 
practice to consider dimensionless ratios of moments of P (q), which are expected to have a 
fixed point (with respect to N) at the critical point, like the (connected) Binder parameter 

mainzi 

and the skewness 



((q-(q)) 2 ) 



S{h : T) = — 3/2 . (6) 

((<1-(<1)) 2 ) 

Here B(h,T) is defined in such a way that it is zero for a Gaussian distribution and one 
for a two equal weights 5-function distribution, whereas the skewness is a measure of P(q) 
asymmetry, which is non-zero in our case of non-zero magnetic field. 

Further evidence for the presence of a phase transition should come from the behavior of 
the spin glass susceptibility 

XsG (^,iv)ociv(F)-¥ 2 ), (7) 

which is expected to diverge, in the thermodynamic limit, when entering in the spin glass 
phase, like 

Xso(h, oo) oc — h -> h\ T , T fixed. (8) 

h - h AT 

According to usual finite size scaling arguments, the finite size behavior is (for h near to ^at) : 
X sg{K N) = N^xsg [N 1/3 (h - h AT )] , T fixed. (9) 

The finite size behavior of the spin glass susceptibility in the SK model on the AT line was 
numerically studied in [IK] , and the scaling « N 1 ^ 3 was checked. 

The non-self-averageness of P(q) is among the many fascinating features of the SK model. 
For instance, considering four replicas of the system, one finds jHIj that Pj{qi2}Pj{qz^) ^ 
P(Qi2)P(qu), whereas the following relations hold: 

P(gi2,g 3 4) = Pj{qi2)Pj{qu) = \p{qi2)P{q^) + \p{qi2)5{qi2 - g 34 ) (10) 

P(qi2,qi3) = PAqi2)Pj(qi 3 ) = \p(qi2)P{qn) + ^P(gi 2 )<%i 2 - q 13 ). (11) 

These kind of relations, which are non-trivially verified if the replica symmetry is broken 
{P(q) $(q ~ qEA))i were recently derived under very general conditions [S^-jSS], such as 
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stochastic stability |2E|- Infinitely many sum rules follow. We consider in particular the 
relations: 



Kmt(h,T) = (q 12 qu)-^(q) ~ ^(q 2 ) = (12) 



(q ) - (q) xsG{n, i ) 6 



1-TT2 1- 



Rl2ls(h,T) = (g 12 g 13 ) - -<g) - -<g2) = Q (14) 



-2 



^ 2 34(^ t) = (SU-^F) 2 -^F) = o (16) 

l^r2 1- 



^m 3 (^) = (g 1 2 2 g 1 2 3 )-^('? 2 r-^g 4 ) = (17) 

Relations R\234:(h,T) and i?i 2 fe 13 (/i, T), which have to our knowledge never been previously 
investigated numerically, are expected to be verified only at non-zero magnetic field, since 
these relations are not invariant under a global flip, and an infinitesimal magnetic field was 
implicit in the derivation of equations (jlOHllj) . On the other hand, relations -R 2 234 (/i, T) and 
i? 2 213 (/i, T) are valid for zero magnetic field also, and were already studied (for 3d and Ad 
Ising spin glasses at zero magnetic field EHIED] ). 

We measured also the following ratios of moments which are non-zero when the system is 
non-self-averaging, and have been introduced for locating the transition point 



G(h,T) 



(q 2 ) 2 


-(q 2 ) 


<<? 4 >- 


- <sT 


((?- 


(q)) 2 ) 2 


<(<?- 


■ <<?» 4 > ■ 


(q 2 ) 2 


- (q 2 ) 2 



G c (h,T) = ^" 1 ^ w/y ' 2 (19) 



A(h, T) = w 1 l q 1 (20) 
(g 2 ) 



A c (h, T) = <fa (fa T to)P>_ (21) 

«9-<9» 2 > 

In the infinite volume limit, G(h,T) is expected to take the constant value 1/3 in the 
glassy phase, because of relation (fTUJ). whereas G c (h,T), A(h,T) and A c (h, T) are non-trivial 
function of h and T, that are zero in the whole paramagnetic phase. 
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These parameters should be particularly useful when the Binder parameter behaves non 
monotonically, taking both positive and negative values, as it is found to happen when there is 
no time- reversal symmetry in the Hamiltonian El EHl El H2] (like in our case of non-zero 
magnetic field) and in systems where the mean field solution is one step replica symmetry 
breaking (121 • Their behavior has been extensively studied and they have been applied to a 
number of models [2E1 12H1 III 1121 EHl III • 

The study of A and A c allows to check if the numerator in G and G c is really non zero 
or if it is approaching zero for increasing volumes together with (or more slowly than) the 
denominator [2E1 EH] • This is obviously not expected to happen in the glassy phase of the SK 
model. As a matter of fact, it was shown j2H] that G and G c should take the universal values 
1/3 and 13/31 respectively at zero temperature for any finite volume Ising system under quite 
general hypothesis, i.e. even if the order parameter is self-averaging. 

The connected parameter G c should be the most effective quantity to look at for locating 
the transition point, since it seems to be the one less affected by finite size corrections to 
scaling j2H]. 

Relations -R1234, #1213 should behave as G and are in principle good candidates for ob- 
taining evidence of the transition. Nevertheless we will show that their finite size behavior, 
as well as the one of G itself, is definitely different from that of G c (at least for the considered 
N values) and that they do not help to locate the transition point. 

The main source of finite size effects is the global reversal of all spins. It does not occur in 
the thermodynamic limit, but when h is exactly zero. It does occur however in a finite volume, 
and as a consequence P(q) develops a tail in the q < region. This tail is significant [To! 12*7] 
even for a size as large as iV = 1024 and a magnetic field value h = 0.3 (at temperature 
T = 0.6). This was observed also in finite dimensional systems fH\ IT*"1 \W\ |2"U] and it is 
expected to strongly affect the scaling of different quantities. In order to reduce its importance, 
one may use [HI [1*5] the "absolute" spin glass susceptibility defined as 



More generally, one can define "absolute" variants of all quantities defined in Equations 
(|Hlffi|). ()12ll2~T]l . In the following we will systematically study the differences between usual 
quantities and "absolute" ones (that will be labeled by the superscript abs), trying to clarify 
which are the most appropriate to look at in order to get evidence for the transition. 

3 Simulations 

We fixed the temperature at T = 0.6, where the AT line corresponds to the critical value 
hAr{T = 0.6) ~ 0.382 We use the magnetic Parallel Tempering algorithm (h-PT), 




(22) 
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described in details in our previous paper [27]. We consider n = 49 replicas, each one at a 
different magnetic field h from a set of n different values both within and without the AT 
lines, from h m i n = —0.6 to h max = 0.6 at equally spaced intervals of Sh = 0.025. Exchange of 
h values between nearest neighbor replicas are allowed with the usual Monte Carlo acceptance 
probability. Moreover the sign of all spins of a replica at h — can be reversed with probability 
1/2. This makes easier the passage from negative to positive h values and vice-versa. 

This h-PT algorithm is ideally suited to obtain the behavior of quantities as function of 
the magnetic field at fixed temperature. However it was found [2JJ that its efficiency rapidly 
decreases when simulating large systems, most probably because of chaos with magnetic field. 
At variance with the case of temperature chaos, the effect of chaos with magnetic field becomes 
evident already on a size of order N = 1024 and Sh ~ 0.15. This means that the phase spaces 
explored by the system at equilibrium at h and h + Sh become quite different when Sh ~ 0.15. 
Therefore N = 1024 is the largest size we could thermalize with this method. 

We perform 50000+50000, 100000+100000 and 300000+300000 h-PT sweeps for N = 64, 
256 and 1024 respectively, the first half of each run being discarded from the statistics. Ther- 
malization was checked by comparing the data obtained in the second part with the ones 
of the second quarter. We simulated four sets of replicas evolving contemporaneously and 
independently (i.e. 49 x 4 = 196 replicas). Data are averaged over 256 disorder configura- 
tions for each system size. Statistical errors are evaluated from (disorder) sample-to-sample 
fluctuations by using the Jack-knife method. 

Looking for evidences of the shape of P(q) with h ^ predicted by the solution of the 
model (Equations ©), we performed a large scale simulation for a system of N = 3200 spins, 
with the usual PT in temperature algorithm, taking n = 38 equally spaced temperature values 
between T min = 0.4 and T max = 1.325, at magnetic field h = 0.3. Results were averaged over 
128 different disorder realization and we performed 400000 + 400000 PT steps for each sample, 
checking thermalization by comparing the P(q) obtained in the second half and in the second 
quarter of the run. 

4 Results and Discussion 

4.1 Energy, magnetization and magnetic susceptibility 

We plot in [Fig. 1] the energy density and the magnetization as a function of h for the different 
sizes considered. In the same figure, we also present data on the mean overlap (q) and on the 
mean absolute value of the overlap (\q\). It is evident from these data that the two quantities 
definitely differ for h as large as h ~ 0.4 (i.e. larger than /zat) for N = 64 and up to h ~ 0.2 
for N = 1024. 
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Figure 1: The energy density (top left), the magnetization density (top right) and the mean 
value of the overlap : (q) and (bottom) at T = 0.6, as a function of the magnetic field, 
for the different system sizes. 
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Figure 2: On the left: usual magnetic susceptibility at T = 0.6, as a function of the magnetic 
field, for the different values of the system size, compared with the predicted infinite volume 
analytic behavior (see text). On the right: the magnetic susceptibility from the absolute value 
of the magnetization. In both plots the two estimates of the susceptibility (|2*3j) and ({231 are 
plotted. 



In [Fig. 2] we plot the magnetic susceptibility computed as 

X(KN) = % (23) 

and as 

x(h,N) = N(3((m 2 ) - (m) 2 ). (24) 

Excellent agreement 1 is observed between the two estimates, showing that good sampling 
has been achieved. We also plot the analytical estimate jjjj (obtained in the so-called Parisi 
approximation for the Landau potential, and for T — > T~): 

X (h,oo) = l-(3/4)^h 4 / 3 , (25) 

which is in quite good agreement with the data for h not too small 2 . In the thermodynamic 
limit, x — /3(1 — (<?)) = ft/P — 1 in the whole spin glass phase for h — > 0. On a finite system 
however one must take into account the symmetry with respect to the flip of all the spins and 
x(h = 0, oo) = f3(l — (q)) = (3 ~ 1.666. Our data for x show clearly the crossover between 

Although formula [251 and 121 are identical mathematically, their Monte Carlo estimates can disagree, if 
the sampling is bad, or if 5h (we use finite difference to estimate I23fl is too large. 

2 A different formula appears in 02] and e.g. ^HJi with the order /i 4 / 3 term multiplied by 7/3. 
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these h = and h 7^ regimes. Figure 2 shows that finite size effects are positive for low h, 
and negative for larger h. This change of sign must happen since the susceptibility has sizable 
finite size corrections whereas the magnetization is (trivially) size independent for h = 0, and 
very weakly size dependent for large h. 

In order to better clarify the situation, we consider also the "absolute" magnetic suscepti- 
bility x abs (see the same figure on the right), obtained from the probability distribution of the 
absolute value of the magnetization. Note that for the "absolute" susceptibility the estimates 
of equations and do not coincide, but when finite size effects are negligible. 

4.2 On the critical behavior at the AT line 

Let us now look at the dimensionless ratios of P(q) moments which should intersect at }iat{T) 
giving evidences for the phase transition (Here Hat — 0.382). We present in [Fig. 3] the 
behavior of the Binder parameter and of the skewness. The original B and S of Equation 
and © on the left of the figure, the "absolute" variants on the right. We see that in all cases it 
would be a hard task to get a clear unambiguous determination of the critical point from the 
data. At variance with the behavior of the h = Binder parameter which is always positive 
and increases continuously for decreasing T's, all four quantities in [Fig 3] (the absolute value 
does not help) display a non-monotonic behavior, taking negative values on a large part of 
the interval. For B and S, the intersection between N = 256 and N = 1024 happens at a 
definitely too small h ~ 0.2, showing the presence of strong finite size scaling corrections (In 
all cases the value A = 64 turns out to be too small to give interesting results). The effect of 
the absolute value is strong, particularly for iV = 64 and 256. It goes in the right direction 
since comparing data for A = 64, 256 and 1024, we now find a monotonous crossing point 
behavior which approach Hat from above, with the intersection between iV = 256 and 1024 
data happening (within the error) at the correct value Hat, the agreement being best in the 
case of the skewness. 

Next we consider in [Fig. 4] the parameters based on P(q) glassy phase non-self- aver ageness. 
The first observation is that in order to obtain some information one has to look at connected 
quantities, since G and A have a definitely different behavior from the others and do not 
cross as one could expect. More in detail, curves for G corresponding to iV = 64 and 256 are 
nearly constant (within the errors) in the whole considered /i-range. For A = 1024, G is in 
agreement within errors with the thermodynamic limit value 1/3 in the glassy phase, and is 
smaller in the paramagnetic phase. For h as large as 0.6 one still finds G ~ 0.28. In the case 
of the EA model at fixed magnetic field the parameter G was found [2H| less appropriate than 
G c . It decreases more slowly when entering in the paramagnetic region and curves for different 
sizes cross at definitely too large temperatures (here quite small sizes with A between 5 and 
64 were considered). 
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Figure 3: The behavior of the Binder parameter B(h, T) (up, left) and of the skewness S(h, T) 
(bottom, left) at T = 0.6, as a function of the magnetic field for the different system sizes. 
Here Hat — 0.382. On the right are plotted the corresponding quantities B abs (h,T) (up, 
right) and S abs (h,T) (bottom, right) obtained from the distribution of the absolute values of 
the overlap. 
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Curves for A are monotonic and approach zero (though with strong finite size corrections) 
at large h, making evident that P(q) is self-averaging in the paramagnetic phase, but they 
do not cross correctly and in particular at h — one finds A(N = 1024) ~ A(N = 256) < 
A(N = 64). 

The connected parameters display a more interesting behavior and they do indeed cross. 
We find a qualitative similar behavior near the transition point for G c and A c . Our statis- 
tical errors are quite large and do not allow a precise determination of the crossing points, 
nevertheless curves for the two smaller sizes seem to roughly intersect at the correct value 
h ~ 0.4 ~ fiAT, whereas data for N = 256 and 1024 appear to cross at a lower /i-value, 
this being more evident in the case of A c . The behavior seems still better when looking at 
the corresponding "absolute" parameters. In particular, there is no irregular behavior for 
h — ► and the statistical errors are definitely smaller, allowing a more precise evaluation of 
the crossing points. Curves corresponding to N = 64 and 256 intersect at h ~ 0.6 whereas 
N = 256 and 1024 are clearly crossing at the right value, h ~ 0.4 ~ Hat- 

To conclude the discussion on these parameters, we note that although G c and A c are cer- 
tainly interesting to look at for getting evidences of the transition (and of non-self-averaging), 
one should not overlook their large statistical fluctuations, much larger than the fluctuations 
of the usual Binder parameter or the skewness, therefore a larger number of samples would 
be required in order to obtain precise measures. 

An even better evidence for the presence of a phase transition comes from the behavior 
of the spin glass susceptibility which is clearly diverging (see [Fig. 5]) when entering in the 
spin glass phase. We see that the behavior of x a sc (plotted in the same [Fig. 5]) is definitely 
different for the largest size considered N = 1024 also. Here the susceptibility seems to 
approach a constant in the small field region, i.e. for h < h m i n (N). Nevertheless h m i n (N) 
is clearly approaching zero for increasing sizes and also in this case one gets evidence for a 
diverging spin glass susceptibility. 

The importance of considering both the usual xsg an d the corresponding "absolute" 
quantity Xsg becomes evident when looking at the scaling plot presented in [Fig. 5]. We find 
strong corrections to scaling in the glassy phase, such that it would be hard to evaluate the 
correct critical point and exponents from these data. On the other hand, corrections to xsg 
are in the opposite direction than those on x a s6i an d therefore it is useful to look at both 
quantities for understanding the true scaling behavior. We also note that data for xsg an d 
Xsq are nearly coincident in the whole relevant interval (i.e. down to h ~ 0.3) for N = 1024, 
which means that when looking at sizes of this order or larger one can expect to find not too 
important corrections to scaling. 
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Figure 4: The behavior of the parameters G(h, T) and A(h, T) (top), G c (h, T) and Gf s {h, T) 
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4.3 The sum rules 



Our next aim is to investigate the validity of the stochastic stability sum rules (jlOH17|) . Some 
are supposed to be non-trivially valid for h < h A r, and trivially valid for h > h A r, namely 
#12345 #12135 #1234 an d #1213- n the other hand relations #J 2 34 an d #12135 that are ratio of 
moments, are supposed to be only valid below the AT line. We will also consider relations 
Ri2U and #1234, obtained from the probability distribution of the absolute values of the 
overlap. 

We first look at the behavior, as function of the magnetic field, of the different terms 
entering the sum- rules, considering the case of N = 1024, where finite size corrections are less 
important. We see from [Fig. 6] that 



-2 



(<?i2<?34) ^ (\q\) 2 ^ (312013) - (?) - (<? 2 > (26) 

for h > 0.4, which is quite reasonable since this is the replica symmetric region where P(q) is 
self-averaging. The tail of P(q) in the negative overlap region has practically disappeared for 
N = 1024, so that (q) ~ (\q\), and we have moreover (assuming a Gaussian distribution with 
variance a 2 ): 

W) = cr 2 + q 2 EA ^q 2 EA =Jq) 2 . (27) 

On the other hand, these quantities are different when entering into the glassy phase, 
though the differences appear small. 

We plot in the same [Fig. 6] relations #1234, R1234 an d #1213- I n order to show that they 
are non-trivially verified we plot also, following the relations 

T 12U (h,T) = Jq-^-^) 2 -^) = (28) 



1- 



T^kiKT) = (|g|) 2 --(|g|) - -<g 2 > = (29) 

" 2^2 



T 1213 (h,T) = ( qi2 q 13 ) - -( q ) — ~{q 2 ) — 0. (30) 

which should also be verified if -R1234, R1234 an d #1213 where trivial. They are clearly not 
verified in the spin-glass phase and accordingly -R1234 an d #1213 are non-trivial in this phase. 
We moreover note that to look at the probability distribution of the absolute value overlap 
is very useful also in this case, since as we already pointed out these relations are derived 
assuming an infinitesimal magnetic field which breaks the global symmetry for inversion of 
all the spins. As a matter of fact, relations -R1234 an d #1213 are no m ore verified as soon as 
h ~ 0.15, where the tail of P (q) in the negative overlap region becomes important also for 
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right) as a function of the magnetic field for 



iV = 1024 (we do not present data for h < 0.1). On the other hand, relation -R1234 appears 
very well satisfied within the errors down to h = 0. 

In [Fig. 7] we present our data for relations -R1234, #1234* an d -R1213 with a finer vertical 
scale than in [Fig. 6], together with data for -RJ234, -R1234 an d #1213- The situation is quite 
clear: for N = 1024 relations -R1234, #1234; #1213 are no ^ verified above the AT line, whereas 
they are satisfied within the statistical error below (up to crossover effects for small h's for 
the non "absolute" quantities) . 

In any event, the change of behavior in the sum rules when going from the h < 0.4 region 
to the other side of the AT line is small. This is trivial in the case of relations -R1234, R123A an d 
#1213 srnce a U terms become very similar (see [Fig. 6]) as we already discussed. This can be 
understood jSO] for R b 123i , R h {$£ and R b 1213 (i.e. the ratios XsJVxsg = xfc^ I Xabs.sc = 1/3 
and XsgVxsg = 1/2 respectively), using the results of jSI]. This R b 's can be calculated from 
the masses Tr, tl and (r^ — r^jn computed in this paper. Rather surprisingly, the ratio 
condition becomes true again in the high field limit, and the R b, s gain back their 1/2 and 1/3 
values. This means that these R b, s have only a very slight variation in the RS phase, with 
probably a minimum, and they are continuous at the AT-line. We note that their behavior is 
very similar to the behavior of G, further confirming that appropriate parameters for getting 
evidence for the transition are the ones which involve connected quantities, such as G c and 
A c . 

In [Fig. 6] we also present the behavior of the different terms entering the relations -R1234 
and #1213- (<7 2 ) 2 i s definitely different from the other terms, and quite surprising remains 
clearly different also on the other side of the AT line, which is to be interpreted as a reminis- 
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cence of non-self-averageness due to finite size effects. 

Finally we plot in [Fig. 8] relations -R1234 an< ^ -^1213- Here finite size effects are less 
important because these relations are valid also in the h — > limit. Nevertheless these 
quantities are compatible with zero within our statistics for h ~ 0.3 only for iV = 1024. 
Also in this case we find only a small difference between the h < 0.4 behavior and the one 
outside the glassy phase. From this point of view it should be recalled that (q 2 ) 2 is definitely 
different from the other terms in the whole h range, which means that data for h > 0.4 are 
far from being in the asymptotic self- averaging regime in which these sum rules should be 
trivially satisfied. As a last remark, we note that from relation -Ri 234 immediately follows the 
expected behavior of the parameter G and that the small differences we observe here between 
the behaviors inside and outside the glassy phase do indeed reflect the fact that G is not an 
appropriate observable to look at for obtaining evidence of the transition. 

4.4 On P(q) for a large size 

Our data for P(q) at h — 0.3 for a system of 3200 spins can be found in [Fig 9] for T = 
0.4, 0.5 and 0.6. The corresponding values of qEA are 0.759, 0.640 and 0.505 respectively. 
We have been very careful in checking that thermalization is achieved for all values of the 
temperature. It is clear from the figure why the asymptotic behavior of this distribution 
has escaped observation up to now. At T = 0.6 only a single peak (corresponding to qEA) 
is visible, with substantial asymmetry (the distribution is wider in the low q side). The 
asymmetry is stronger for T = 0.5, but there is still no sign of the low q peak 3 . Only for 
T = 0.4 does one see the expected continuum on the left of the self-overlap peak, with some 
indication of the low q peak at a location in agreement with the value jlH] qmin — 0.44. It 
should be noted that the peak corresponding to q min is predicted to be broader than the one 
corresponding to qEA • This explains why we do not observe this minimum overlap peak. 

5 Conclusions 

We performed numerical simulations of the SK model in a magnetic field at temperature 
T = 0.6, both in the glassy phase and above the AT line. We used a modified version of the 
PT algorithm in which the system is allowed to move between a chosen set of magnetic field 
values, an algorithm well suited for our purpose. 

We measured quantities such as the magnetic susceptibility, which turns out to be in 
agreement with the predicted analytical behavior of [UJ as function of h. 

3 This result disagree with [T^l, where a low q peak is found, using the Metropolis algorithm with 100000 
sweeps for equilibrium and only 20 disorder samples. 
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Figure 9: The behavior of the probability distribution of the overlap P(q) with h = 0.3 and 
temperatures T = 0.4, 0.5 and 0.6 respectively, for the large size N = 3200. 
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Dimensionless ratios of P(q) moments such as the Binder parameter and the skewness 
display a non-monotonic behavior making difficult to get a clear determination of the tran- 
sition point on the AT line. Also ad hoc parameters for locating replica symmetry breaking 
transitions, based on the non-self-averageness of the order parameter, are considered. The 
connected ones turn out to be effective for locating the transition. 

An even better evidence for the transition comes from the divergence of the spin glass 
susceptibility, though its scaling behavior is affected by strong finite size corrections. 

We also investigate the behavior of various quantities defined in term of the probability 
distribution of the absolute value of the overlap. This allows to reduce the finite size effects due 
to the long tail of P (q) in the negative overlap region. As a matter of fact, the dimensionless 
parameters turn out to behave better in this case, the crossing points being nearer to the 
correct critical value. It is interesting to note that the usual and "absolute" susceptibilities 
have corrections of opposite signs. 

Moreover we studied the behavior of some sum rules (related to stochastic stability) in- 
volving overlaps between three and four replicas. We found strong finite size corrections 
particularly for those sum rules that are valid only at non-zero magnetic field, and it turns 
out to be particularly appropriate to look at "absolute" quantities in this case. They are 
satisfied within our statistical accuracy for N = 1024 in the glassy phase. On the other hand, 
they would not be good indicators for the transition, since their behavior change very slightly 
when crossing the AT line, being still nearly verified also for h > Hat, some trivially (all the 
terms become very similar) others non-trivially. 

Finally we presented data for P(q) in magnetic field, which show how slowly the shape 
predicted by the RSB solution develops on a large system. 
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